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Abstract. An effective field theory is derived that describes the low- frequency spin dynamics in 
the low-temperature orthorombic phase of LoiCuOa- Restricted to a single CuOi layer the effective 
theory is a simple generalization of the relativistic nonlinear a model to include all spin interactions 
allowed by symmetry. Incorporating a weak interlayer interaction leads to two coupled nonlinear a 
models which provide an efficient description of the complete bilayer dynamics. Particular attention 
is paid to the weak-ferromagnetic and spin-flop transitions induced by external magnetic fields. The 
main features of the observed (covert) weak ferromagnetism are thus accounted for in a straight- 
forward manner but some of the finer theoretical predictions would require further experimental 
investigation. The derived framework is also suitable for the study of the structure and dynamics 
of magnetic domains in undoped Lo^CuOa- 

PACS. 75.10.-b General theory and models of magnetic ordering — 74.72.Dn. La - based cuprates 



I. INTRODUCTION 

The magnetic properties of La-iCuO^ have been ex- 
tensively studied during the last decade Q). In a 
crude approximation this system is described by a two- 
dimensional (2D) isotropic Heisenberg antiferromagnet. 
However the orthorombic distortion of the crystal that 
takes place below 530i^ induces anisotropic spin cou- 
plings, the most important of which is a Dzyaloshinskii- 
Moriya (DM) anisotropy that should lead to spin 
canting and thus weak ferromagnetism. Yet a small anti- 
ferromagnetic interlayer coupling forces successive Cu02 
layers to cant in opposite direction and the induced 
weak moments average to zero in the absence of external 
fields. Nevertheless the sustained experimental as well 
as theoretical effort of Thio et al. d-[| established that 
La^CuO^ is indeed a covert weak ferromagnet. 

The derived phcnomcnological picture was further 
probed by symmetry analysis and by a fresh and 
challenging look at the DM anisotropy due to Kaplan M 
and Shekhtman, Entin-Wohlman, and Aharony |]Io| , |l if . 
The main outcome of the latter work is sometimes re- 
ferred to as the KSEA anisotropy and is the subject of 
current experimental investigation in the related context 
of helimagnetism |12] . 

Our purpose is to systematize the above developments 
into a simple field theoretical framework that should fa- 
cilitate further work on this interesting subject. The rele- 
vance of effective field theories became apparent through 
standard hydrodynamic approaches |P3|,[h| which eventu- 
ally led to a successful description of the isotropic Heisen- 
berg antiferromagnet in terms of a relativistic nonlinear 
a model Ej^j ). A similar approach has been employed 
for the study of the dynamics of domain walls and re- 
lated topological magnetic solitons in conventional weak 
ferromagnets ]T^ ]. This issue has remained largely un- 
explored in the cuprates, apparently due to the hidden 



nature of weak ferromagnetism and a corresponding lack 
of a complete field theoretical description. 

In Section 2 we repeat the symmetry analysis to ob- 
tain the most general spin Hamiltonian involving nearest- 
neighbor (nn) interactions. Based on this Hamiltonian 
we proceed with the derivation of an effective field the- 
ory valid at low frequencies in two steps. First, in Section 
3, we derive a suitable extension of the 2D nonlinear a 
model which would be appropriate for the description of 
a single CuC>2 layer provided that the interlayer coupling 
could be neglected. Second, in Section 4, we incorporate 
a weak interlayer interaction to obtain the final effective 
theory in the form of two coupled nonlinear er models. 
The main phenomenological implications are also worked 
out in Section 4, including a discussion of the observed 
weak- ferromagnetic Q] and spin- flop Q transitions. Our 
conclusions are summarized in Section 5 by listing those 
issues that seem to deserve further attention. In the Ap- 
pendix we derive the most general next-nearest-neighbor 
(nnn) in-plane interaction compatible with symmetry. 

II. SYMMETRY CONSTRAINTS 

The crystal structure of La^CuO^ has been discussed 
on several occasions and the accumulated information is 
freely used in this section. In particular, we found the 
expositions of references []iq , |l9l lucid and informative. 
The crystal undergoes a structural phase transition from 
a tetragonal (I4/mmm) phase at high temperatures to 
an orthorombic (Bmab) phase below 530-fiT. Throughout 
this paper we confine attention to the low-temperature 
orthorombic (LTO) phase. 

The relevant space group Bmab is usually listed as 
Cmca in standard tables of crystallography using a 
slightly different choice of conventions. Our conventions 
are illustrated in Figure 1 which depicts the orthorombic 
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unit cell stripped of all but the magnetic sites; i.e., the 
positions of the spin s = 1/2 Cu 2+ ions. Dashed lines 
in this figure join nn neighbors within each Cu0 2 plane 
and are the descendants of the original tetragonal axes 
which are no longer orthogonal due to the orthorombic 
distortion. In a first approximation, the magnetic ground 
state is such that spins at sites denoted by A and A point 
along the positive (negative) b axis, while spins at B and 
r point along the negative (positive) b axis |2(|. How- 
ever, when anisotropics and a weak interlayer coupling 
are taken into account, each spin suffers a slight cant- 
ing which leads to four inequivalent magnetic sites and 
a corresponding four-sublattice picture. In the follow- 
ing, spin vectors are denoted by their standard symbol 
S, or by A,B,T, A when a distinction among the four 
sublattices becomes necessary. 
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FIG. 1. The orthorombic unit cell of La2CuOi stripped 
of all but the Cu atoms denoted by solid circles. Dashed 
lines join nn neighbors within each CuOi plane and are the 
descendants of the original tetragonal axes. The four inequiv- 
alent magnetic sites are labeled by A, B, T and A. 

In order to list the symmetry elements of the space 
group we begin with the most general primitive transla- 
tion 



where e a , e b and e c are unit vectors along the crystal 
axes, a = 5.35A, b = 5.40A and c = 13.15A are the lat- 
tice constants, and (a, /?, 7) is a set of integers that may 
also be used to label the relative position of a unit cell 
through the Cartesian coordinates x = aa, y = (3b and 
z = 7c. We further consider the two fractional transla- 
tions 



r = -{ae a + ce c ), 



r = -(ae a 



be h 



(2) 



where r is in itself a symmetry operation. The symmetry 
group of the unit cell is then written symbolically as G = 
Go + txGq where Gq contains the eight elements 



E, I, a a , a' b , a' c , C2a, C' 2b , C' 2 



2c- 



(3) 



(3be b 



■jce c 



(1) 



Here E denotes identity, / inversion about a Cu site, 
a a ,ab,cr c reflections about the planes x = a/2, y = fe/2, 
z = c/2, and C 2a , C 2 b, G 2c 180° rotations around the axes 
that emanate from the center of the unit cell as shown in 
Figure 1. Primed elements in equation (3) must be com- 
plemented by the special fractional translation r' which is 
not in itself a symmetry operation. Finally we note that 
all elements of G can be generated from the fundamental 
set (T,I,a a ,C 2c ) by suitable group multiplications. The 
above four elements have been used to derive the spin 
Hamiltonian described in the remainder of this section 
without presenting the detailed symmetry arguments. 

We first discuss a 2D restriction that would be appro- 
priate for the description of a single Cu0 2 layer if the 
interlayer coupling could be neglected. The special sym- 
metry element r of equation (2) maps AB bonds in the 
basal plane to TA bonds in the middle plane. Therefore 
it is sufficient to consider for the moment only the basal 
plane. We assume that the Hamiltonian contains terms 
that are at most quadratic in the spin operators, even 
though quartic terms are occasionally invoked to explain 
some features of multimagnon dynamics |2lf| . Including 
all nn interactions within the plane the Hamiltonian is 
written as the sum of four terms : 

W = W e + Wdm +W A + W Z , (4) 

which we now describe in detail. The first term 

W E = ^2 Ju(S k ■ S t ) (5) 

<kl> 

contains the isotropic exchange interaction over nn in- 
plane bonds denoted by < kl >. Symmetry requires that 

Jki = J, for all in— plane nn bonds. (6) 

The second term 

W DM = Gkl ■ ( Sfc x S ') ( ? ) 

<kl> 
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is the standard (antisymmetric) DM anisotropy [§,1). 
The vectors are severely restricted by symmetry and 
take only two distinct values: 



Vp+i 



Di =De a + D'e b , 



D 2 = De a - D'e b , 



(8) 



which are distributed over the 2D lattice as shown in 
Figure 2 where a sign alternation ±Di and ±D 2 on op- 
posite bonds is also displayed. This alternation together 
with the specific form of DM vectors (8) were derived 
by Coffey, Bedell, and Trugman 0. Actually the DM 
vectors given in the above reference differ from those of 
equation (8) by a 45° rotation, apparently because they 
were referred to the original tetragonal axes. Since the 
latter axes are not exactly orthogonal in the LTO phase, 
statement (8) should be viewed as slightly more precise. 
The third term 



w * = \ E +4$) 



(9) 



<kl> i,j 



encompasses all "symmetric" exchange anisotropics over 
nn in-plane bonds. Again the matrices are restricted 
by symmetry to two possible values: 



Ki = 



K aa 


K ab 





K ab 


K bb 














K aa 


-K ab 






-K ab 


K bb 




(10) 








K cc 





which are distributed as shown in Figure 2. The above 
matrices may be taken to be traceless (K aa + K bb + K cc = 
0) because the isotropic component of the exchange in- 
teraction has already been accounted for by equation (5). 
Finally the fourth term 



W z 



(ii) 



is simply the Zeeman interaction produced by an external 
field H. 

The symmetry analysis was further extended to include 
all nnn in-plane interactions; namely, couplings along the 
diagonals of the Cu plaquettes which are parallel to the 
orthorombic a and b axes. The resulting additions to 
the Hamiltonian are summarized in the Appendix. How- 
ever, although we have kept track of the nnn interactions 
throughout our analysis, the corresponding results will 
not be included in the main text because we wish to 
keep the exposition reasonably simple. But we will oc- 
casionally mention potential contributions from the nnn 
interactions. 

We next turn our attention to possible microscopic 
mechanisms that produce the various anisotropies. A 
good starting point is the KSEA Hamiltonian [fill 
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FIG. 2. Illustration of the distribution of the DM vectors 
±Di and ±D2, and of the symmetric exchange matrices K\ 
and Ki , on a portion of the basal plane. The spins A and B 
in a given dimer are labeled by a pair of indices (a, j3) that 
advance along the orthorombic axes a and b not shown in the 
figure. 



w = y, [ (J™ 

<kl> 
1 



2Jfci 



4J H 

S fc .D fcl )(D w -S I )] > 



)(S fe -S z ) + D fci -(S fc xS,) 

(12) 



obtained by carrying the perturbative treatment of spin- 
orbit interaction to second order. The second-order terms 
in equation (12) are usually neglected by comparison 
to the linear DM anisotropy. However their effect can 
be far more subtle noting that the total contribution 
from each specific mode can be brought to a completely 
isotropic form by rotating the spin operators and S/ 
around the D^; axis with angles —6ki and 8ki such that 
tan6>fcj = \Dki\ 2 /2Ju- One would naively conclude that 
Hamiltonian (12) does not predict spin ordering of any 
kind; in particular, weak ferromagnetism. Nevertheless 
spin ordering may occur as a result of frustration when 
the collective effect from all bonds is taken into account. 
This delicate line of reasoning |ll[] led to some inter- 
esting predictions such as the existence of weak ferro- 
magnetism in the LTO phase of pure LaiCuO^ but its 
absence in the low-temperature tetragonal (LTT, space 
group PAiilncm) phase that may be realized in, say, Ba- 
doped La2~ x Ba x CuO,±. 

It should be expected that the exchange anisotropy ab- 
stracted from equation (12) is consistent with the sym- 
metry statement (10). Indeed, taking into account that 
Jki = J for all in-plane bonds and the specific form of 
the DM vectors from equation (8), the (traceless) KSEA 
anisotropy is found to be a special case of (10) with 



Kaa = 



K c ,= 



2D 2 - D r 



6J 



D 2 + £>'■ 



6J 



K bb = 



K ab = 



2D r 



D 2 



6 J 
DD' 



2 J 



(13) 
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Of course, spin-orbit interaction is not the only source of 
anisotropy and, in fact, Coulomb exchange produces an 
Ising-like term Jll],p2| described by 



K aa = X -K = K hh , Koc = ~K, (14) 

which should be added to the corresponding elements of 
equation (13). Therefore we proceed with caution using 
the most general anisotropy given by equation (10) and 
return to further discussion of the special cases (13) and 
(14) after the dust settles, in Section 3. 

The remainder of this section is devoted to a brief ex- 
planation of our main strategy. At this point we invoke 
the classical approximation, the validity of which will not 
be questioned until Section 4. The 2D dynamics will thus 
be described by the Landau-Lifshitz equation p3| in the 
form given by Gilbert p4| : 



dA 
~dt 



7(A x 



dA , 

8t ' 

9B, 



= AxF, 



A = 



dW 
9A : 



(15) 



<9B ,„ <9B . „ dW 

- +7 (Bx-)=BxF B , F B = - — , 

where 7 is a dissipation constant. The spin variables A 
and B are treated as classical vectors of length s, and the 
lattice indices (a, (3) displayed in Figure 2 are suppressed 
in equation (15) for notational convenience. 

The suppression of indices is more than symbolic if one 
wishes to study only the homogeneous spin dynamics in 
the presence of a spatially uniform field H. Then spins 
assume only two distinct values A and B, one for each 
sublattice, and the effective fields Fa and Fb in equation 
(15) may be derived from the much simpler Hamiltonian 

W c = 4[J(A • B) + D(A b B c - A c B b ) 

+K aa A a B a + K bb A b B b + K CC A C B C ] 
-H-(A + B), (16) 

where (A a , A b , A c ) and (B a ,B b ,B c ) are the Cartesian 
components of the spin vectors A and B along the or- 
thorombic axes. A notable fact is that neither D' nor 
K ab enter equation (16). If we further restrict the diag- 
onal anisotropies to the Ising form (14), equation (16) 
yields the unit-cell Hamiltonian employed by Thio et al. 
to account for a wide range of experiments The 
sign alternation of the DM vectors on opposite bonds |t]] 
is crucial for the description in terms of a unit-cell Hamil- 
tonian and was thus implicitly assumed in reference [Q. 
Lack of sign alternation would lead to spiral magnetic 
order or helimagnetism JT^ . 

The information accumulated so far is employed in Sec- 
tion 3 to study the 2D dynamics of a single layer. One 
should recall that the dynamics of the middle plane is 
completely isomorphic and may be obtained by the sim- 
ple substitution (A,B) — > (T, A). The complete 3D dy- 
namics including interlayer interactions will be studied 
in Section 4. 



III. DYNAMICS OF A SINGLE LAYER 

The unit-cell Hamiltonian (16) may be analyzed 
through the classical Landau-Lifshitz equations (15) to 
furnish explicit predictions for the characteristic magnon 
frequencies and the corresponding dynamic susceptibili- 
ties, following the early treatment of orthoferrites by Her- 
rmann [^,^6|. Considerable simplifications are invoked 
along the way by appealing to the phenomenological fact 
that anisotropies and the applied field are much smaller 
than the isotropic exchange constant: 



d,k,h<s:j. 



(17) 



Simply stated our task is to derive an effective low- 
frequency dynamics in which the strong inequalities (17) 
are taken into account from the outset. And, at little ex- 
tra cost, we aim to go beyond the homogeneous dynamics 
to include spatial variations within a complete continuum 
field theory in the form of a nonlinear a model. 

One must thus deal with the complete discrete Landau- 
Lifshitz equations (15) where A = A Q _p and B = ~B a ,[3 
are the spin vectors in a generic dimer labeled by two 
integers (a, (3) that advance along the orthorombic a and 
axes, as shown in Figure 2. The corresponding effective 
fields and F b are then constructed by closely follow- 
ing the prescription for the complete microscopic Hamil- 
tonian (4) given in Section 2. The resulting algebraic 
expressions are somewhat lengthy but the actual deriva- 
tion is a straightforward adaptation of a direct method 
employed earlier for the study of an ID model weak fer- 
romagnet [^7j and a 2D antiferromagnet [^8|. We thus 
suppress the algebraic details and state the final result 
which can be explained in a simple manner. 

Rationalized space-time variables are defined from 



V 



aS, 



2s8Jt, 



(18) 



where 8 is a dimcnsionless scale whose significance will 
be made precise as the argument progresses. The actual 
distances along the a and b axes are then given by x = 
r]a/5 and y — £b/5 while frequency is measured in units 
of 2sS J. We also define rescaled parameters grouped into 
two categories. The DM anisotropy and the applied field 
are scaled linearly with 5 : 



d = 



2D 



H 



(19) 



SJ' 2s8J' 
whereas diagonal anisotropies are scaled quadratically: 



SKaa 
8 2 J 



Kb 



SK bb 
S 2 J ' 



&Kcc 
S 2 J 



(20) 



Note that the parameters D' and K ab do not appear in 
the above list because they eventually drop out of the 
effective low-frequency dynamics, for essentially the same 
reason they do not appear in the unit-cell Hamiltonian 
(16). 
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Concerning the field variables a transparent formula- 
tion is obtained in terms of the "magnetization" m and 
the "staggered magnetization" n which are defined by 



Kb 



K b 



(29) 



1 



(A + B), 



1 



(A-B), 



(21) 



and satisfy the classical constraints m ■ n = and m 2 + 
n 2 = 1. The basis for the derivation of an effective field 
theory is that the strong inequalities (17) imply |m| <C 
|n|. Indeed, a consistent low- frequency reduction of the 
Landau-Lifshitz equation is obtained by treating m as a 
quantity of order S, whereas the staggered moment n and 
the rescaled variables (19)-(21) are assumed to be of order 
unity. Then, to leading order, the classical constraints 
reduce to 



0. 



1. 



(22) 



and m may be expressed in terms of n through the alge- 
braic relation 

m=-[nx (n + d- nxh)-(n, + n^)], (23) 

where d = de a , the dot denotes differentiation with re- 
spect to r, and subscripts differentiation with respect to 
the spatial coordinates r\ and £. The various pieces of 
equation (23) have appeared in earlier studies of conven- 
tional weak ferromagnets 17 27j] and antiferromagnets 
p8f and will not be discussed here in detail. 

The main point is that the dynamical equations may 
now be stated entirely in terms of the staggered moment 
n which satisfies the nonlinear a model 



nx(f + Aft) = 0, n 2 = 1, 



(24) 



where A = is a rescaled dissipation constant, and 

the effective field f may be derived from an action prin- 
ciple : 

5A f 
f = — — , A = / Ld v d£dT, (25) 

on J 

where A is the action and L the corresponding La- 
grangian density 



(26) 



L = L Q — V. 
Here Lq is the "free" Lagrangian 

L = i(h 2 -n 2 -n 2 ) + b,(nxri) (27) 

and V is the "potential" 

V = (h x d) • n + ~(n • h) 2 + l -{a\n 2 a + a\n% (28) 

where (n a , rib, n c) are the Cartesian components of n and 



are the final effective anisotropy constants. 

If the applied field were absent (h = 0) the derived 
field theory would be relativistically invariant. The "ve- 
locity of light" is equal to the magnon velocity in the 
corresponding isotropic antiferromagnet and is scaled out 
of equation (27) thanks to the use of rationalized units. 
Recalling that x = rja/S, y = £b/5, and r = 2sdJ, the 
actual magnon velocities along the a and 6 directions 
are V a = 2saJ and Vb = 2sbJ. The predicted slight 
anisotropy Vb/V a = 1.01 cannot be discerned in current 
experiments which yield an average spinwave velocity 
V sw ~ V a ~ Vb « 850meVA for pure La-iCuO^ samples 
Q. Using an average lattice constant a sw b — 5.375A one 
may extract a classical value for the exchange constant 
J = 158meV which differs from the usually accepted 
,/ = ISbmeV by the calculated quantum-renormalization 
factor 1.18. 

We are now in a position to make contact with the var- 
ious special limits of the exchange anisotropy discussed 
in Section 2. The pure Ising anisotropy (14) leads to 
di = d and a 2 = 8K/S 2 J which correspond to the min- 
imal model of Thio et al. In principle, the KSEA 
anisotropy (13) could lead to a\ ^ d. However an es- 
timate of the DM parameters D and D' within a tight- 
binding model (ll) yields the near equality \D\ w \D'\ 
which simplifies the diagonal elements of equation (13) 
to an Ising-like form, and thus a\ ~ d, whereas the off- 
diagonal element K a b does not appear in the effective 
theory. Another mechanism for a\ ^ d is provided by the 
diagonal nnn exchange anisotropies (see the Appendix) 
but we do not have a way to theoretically estimate those 
parameters. In short, departures from the strict equality 
di = d should be expected but the near equality a\ ~ d 
is likely a good assumption. In any case, our theoretical 
calculations will be carried out in terms of the three pa- 
rameters (ai,a,2,d) and further discussion of this issue is 
postponed to Section 4.3. 

Next we discuss the field-dependent terms in equa- 
tions (27) and (28). First we note the well-known fact 
that the last term in the free Lagrangian (27) breaks 
Lorentz invariance |Tt|| . Incidentally, a mild breakdown 
of Lorentz invariance is also induced by the off-diagonal 
nnn exchange anisotropy discussed in the Appendix. The 
potential (28) contains two distinct contributions from 
the external field. The term (n • h) 2 is an easy-plane 
anisotropy due to the tendency of the two spins in a given 
dimer to antialign in a direction nearly perpendicular to 
the applied field. More subtle is the term (h x d) ■ n 
which couples the external field to the antisymmetric DM 
anisotropy. The existence of such a term was anticipated 
by Andreev and Marchenko p9[ in their phenomenolog- 
ical treatment of conventional weak ferromagnets based 
on symmetry. Although this term is often dismissed in 
the literature [0, it was shown to be important for a 
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proper understanding of domain- wall dynamics Effl. In 
fact, the present work will provide ample evidence for the 
crucial importance of such a term in every aspect of weak 
ferromagnetism. 

Whereas the complete physical picture cannot be es- 
tablished until we incorporate the interlayer coupling, in 
Section 4, the remainder of this section is devoted to the 
derivation of some basic consequences of the single-layer 
dynamics. Applications carried out in this paper pertain 
to homogeneous spin dynamics. One may then neglect 
spatial gradients to write 



m = - [n x (h + d — n x h)] , 
L = ifi 2 + h-(nxh), 

V = (h x d) • n + i(n • h) 2 + \{a\nl + a 2 2 n 2 c ), (30) 

which describe the low-frequency dynamics associated 
with the unit-cell Hamiltonian (16). 

For our current demonstration we assume that the field 
points along the c axis (h a = = hb, h c = h) and thus 
the potential reduces to 



V = hdn b + - [afnl + (aj + h 2 )n 



2 1 

cJ ! 



(31) 



whose local minima are given by n = (0, =pl, 0) and the 
corresponding magnetization is computed from the first 
equation in (30) applied for static n; i.e , n = 0. Hence 
the two ground-state configurations are described by 



n = Te b , 



m = —{h ± d)e c 



(32) 



and are degenerate at zero field. For h > the upper 
sign yields the absolute ground state, whereas the lower 
sign corresponds to a metastable local minimum with 
higher energy. Equation (32) makes it explicit that a 
weak ferromagnetic moment develops along the c axis 
even in the absence of an applied field. 

The computation of the spectrum of normal frequen- 
cies is now straightforward. In terms of the standard 
angular variables 

n a + irib = sin O e J * , n c = cos0, (33) 

the free Lagrangian is written as 

L = i(9 2 +sin 2 9 $ 2 ) + /isin 2 © $, (34) 

and a corresponding angular parametrization of the po- 
tential is obtained simply by inserting equation (33) in 
(31). 

The ground-state configurations are then given by = 
■|, $ = qFf ) while small fluctuations are accounted for 
by making the replacements O = j — 9, $ = =f(^ — <p) 



and keeping terms that are at most quadratic in 9 and 
<f>. Also omitting a trivial additive constant and total 
derivatives one finds that 



L = Ln 



-i(c 2 ±( /) 2 +c 2 ± 2 ), (35) 



where 



)\± — a\ ± hd, 



uj 9 ± = a 2 ± hd + h 



(36) 



are the (squared) characteristic magnon frequencies for 
in-plane and out-of-plane fluctuations, respectively. 

The actual experimental values of the calculated 
"magnon gaps" will be discussed in Section 4. But first 
we must fill a gap in our theoretical arguments concerning 
the choice of the scale parameter S in equations (19)-(21). 
Although S plays an important role in ascertaining the 
relative significance of the various terms that arise dur- 
ing the low-frequency reduction of the Landau-Lifshitz 
equations, consistency requires that all physical predic- 
tions be independent of 5. The magnon gaps (36) pro- 
vide a good illustration of this point by recalling that 
the unit of frequency is 2s6 J. When the right-hand sides 
of equations (36) are multiplied by (2sdJ) 2 the resulting 
expressions are indeed independent of 5 and contain only 
the original microscopic parameters in suitable combina- 
tions. But one may exploit 6 to choose more convenient 
rationalized units as discussed in Section 4.3. 

This section is completed by quoting an explicit re- 
sult for the corresponding dynamic susceptibilities. The 
method of calculation proceeds roughly as follows. Sup- 
pose that the system is found in its ground state in the 
presence of a constant field h along the c axis. A weak AC 
field of arbitrary direction and frequency 10 is then turned 
on and the system eventually relaxes into a steady state 
oscillating with the same frequency thanks to the dissi- 
pative term present in equation (24). Once the steady 
state for the staggered moment n is determined from a 
perturbative solution of the nonlinear a model, in the 
weak-AC-field limit, the magnetization m is calculated 
from the first equation in (30) to yield the sought after 
susceptibilities: 



Xaa — . 



Xbb 



1 



iXui 



{d + hf 



4o, 2 



i\uj 



Xcc - 4 



iXui 



Xab 



Xa 



~Xba 



iS (d + h)u> 



4 a; 2 



iXuj ' 



Xc 



Xbc 



Xcb = o, 



(37) 



where u>x and L02 are the magnon gaps (36) restricted 
to the upper (+) sign; i.e., the gaps above the absolute 



6 



ground state when h > 0. The susceptibilities for the 
metastable ground state can be obtained by the universal 
replacement h — > —h. Since the two types of ground state 
mix in the presence of the interlayer coupling, the results 
of equations (36) and (37) should be interpreted with 
caution. 



IV. INTERLAYER COUPLING 

Although interlayer interactions are expected to be 
weak, they are important for a proper understanding of 
spin dynamics in La2CuO^. We have thus extended the 
symmetry analysis of Section 2 to include nn interlayer 
couplings on bonds that are parallel to either the ac or 
the be plane. For each Cu atom there exist eight out-of- 
plane neighbors, four in the plane above and another four 
in the plane below. Symmetry requires that the isotropic 
exchange interaction is described by an exchange con- 
stant Ji for bonds that are parallel to the ac plane, and 
a second exchange constant J 2 for bonds parallel to the be 
plane. These two constants would be equal in the tetrag- 
onal (IA/mmm) phase but are different in the LTO phase 
due to the orthorombic distortion. 

The exchange constants are expected to individually 
satisfy the strong inequalities 



Ji,Jz < J, 



(38) 



in view of the fact that the length of out-of-plane nn 
bonds is significantly larger than the length of in-plane 
nn bonds. It is thus reasonable to assume that out- 
of-plane anisotropies (symmetric or antisymmetric) can 
be safely ignored because they are expected to be even 
weaker. In fact, we have worked out the form of all such 
anisotropies compatible with symmetry to convince our- 
selves that they do not bring in potentially new elements. 

Therefore the three-dimensional (3D) unit-cell Hamil- 
tonian can be written as 



W, 



3D 



W C (A,B) + W C {T,A) 



W mt (A,B,T,A), 

(39) 



where the first term is the 2D unit-cell Hamiltonian (16), 
the second term is obtained by the simple substitution 
(A, B) — > (r, A), and Wi nt contains the isotropic inter- 
layer interactions. Simple inspection of Figure 1 leads 
to 



4Ji (A ■ T + B ■ A) + 4J 2 (A • A + B • r). (40) 



This form is somewhat more involved than the one in- 
voked by Thio et al. Q but will eventually lead to the 
same physical picture. 

We must now reformulate the strategy of Section 3 by 
introducing two pairs of variables 



mx = — (A + B), ni 

Is 

m 2 = — (r + A), n 2 

2s 



B) 



(41) 



which would satisfy two identical copies of the 2D non- 
linear a model derived in Section 3 if the interlayer in- 
teraction (40) were neglected. The latter induces a cou- 
pling between the two copies which is especially simple 
if we invoke the strong inequalities (38). Then the mag- 
netizations mi and m 2 are not directly affected by the 
interlayer coupling; i.e., 



mi = ~[ni x (ni + d-ni x h)], 
m 2 = - [ n 2 x (ri 2 + d - n 2 x h) ] , 



(42) 



but the staggered moments ni and n 2 satisfy a coupled 
dynamics described by the total Lagrangian 



L = L Q - V, 

Lq = Lqi + Lq 2 , 



V = Vi + V 2 + V v . 



(43) 



where Lqi and L02 are two identical copies of the free 
Lagrangian of equation (30) applied for n = ni and n = 
n 2 , respectively ,Vi and V2 are similar copies of the 2D 
potential, and 



V12 = p 2 (ni • n 2 ), 



8(Ji - J 2 ) 
5 2 J ' 



(44) 



is an effective interlayer potential. 

The further inequality J\ > J 2 implied by the nota- 
tion of equation (44) is simply an assumption consistent 
with phenomenology. To illustrate this assumption we 
consider the ground-state configuration(s) at zero exter- 
nal field. In this special case the absolute minimum of 
the total potential V of equation (43) is achieved when 
each term V\, V2 or V\% assumes its least possible value. 
Specifically, ni = — n 2 and 



n 2 = ±e b , 



mi = ±—de c 

4 

S , 

m 2 = T^de c 



(45) 



The spin configuration that corresponds to the upper sign 
is depicted in Figure 3, whereas the second ground state 
is obtained by reversing the sign of all spins and carries 
the same energy. In either case, the average total mag- 
netization m = i (mi + m 2 ) vanishes and thus explains 
the term "covert weak ferromagnet" often employed to 
describe La 2 Cu04. More involved spin configurations 
arise in the presence of external fields and are described 
in Section 4.1 and 4.2. 
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n l 
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FIG. 3. Schematic representation of the ground-state con- 
figuration given by equation (45) with the upper sign. A 
second (degenerate) ground state is obtained by reversing the 
signs of all spins. The canting angle and the corresponding 
magnitude of mi and m.2 are greatly exaggerated for purposes 
of illustration. 

This general description of the effective 3D dynamics 
is completed with three remarks : 

(a) One would think that a more efficient formulation 
could be obtained in terms of the Dzyaloshinskii combi- 
nation of fields H 

m=-(m 1 +m 2 ), n=-(ni+n 2 ), 

P=^(mi-m 2 ), q=^(ni-n 2 ), (46) 

which arise naturally in a proper four-sublattice descrip- 
tion of orthoferrites (26). The latter are conventional 
weak ferromagnets with strong interlayer couplings which 
lead to |m|, |p|, |q| <C |n|. One can then show that the ef- 
fective low-frequency dynamics may be formulated with 
a single order parameter, the total staggered moment n, 
and the remaining auxiliary fields m, p and q are ex- 
pressed in terms of n through algebraic relations similar 
to equation (42). However, such a formulation becomes 
singular when the interlayer exchange constants are much 
weaker than the intralayer one, as is assumed in equation 
(38). As a result, the current formulation in terms of two 
coupled order parameters ni and n 2 is more suitable. 

(b) The assumed strong inequalities (38) also imply 
that gradient terms of any kind produce negligible cor- 
rections to the effective interlayer coupling. Therefore the 
homogeneous 3D dynamics described by equations (42)- 
(44) may be generalized to a complete continuum field 



theory simply by extending the free Lagrangians Lot an d 
L02 to include 2D spatial gradients according to equation 
(27) applied for n = ni and n = n 2 , respectively. The 
resulting field theory is essentially 2D and the only trace 
of 3D dynamics is the bilayer coupling (44). 

(c) The interlayer potential (44) vanishes in the tetrag- 
onal phase (Ji = J 2 ). Hence finer anisotropic couplings 
as well as quantum fluctuations are necessary to explain 
spin ordering pfj| in, say, 5Y 2 Cu0 2 C7 2 which remains 
tetragonal down to temperatures as low as 10 K |3lj]. A 
field theoretical description of these finer couplings is not 
available at present and thus the isotropic 2D nonlinear 
a model is the main tool for investigation of tetragonal 
compounds f32|| . 

A. Weak-ferromagnetic transition 

We return to the problem posed in the concluding 
paragraphs of Section 3 and now address it within its 
proper 3D context. When a field of strength h is applied 
along the c axis the total bilayer potential is given by 

V = p 2 (ni • n 2 ) + hd(n u + n 2 b) 

+ 5 [*?(»?, + O + (4 + h 2 )(n 2 lc + n 2 c )], (47) 

where ni = (m 01 nib, ni c ) and n 2 = (n 2a , n 2b , n 2c ) are 
the two order parameters expressed in Cartesian com- 
ponents. Simple inspection of the potential, taking into 
account that a% < a 2 , suggests that its minima are such 
that ni c = = n 2c . One may then parametrize the re- 
maining components as (nx a ,nib) = (cos $i, sin $i) and 
(«-2a,«2&) = (cos $ 2 , sin<I> 2 ) to obtain 

V = p 2 cos($i - $ 2 ) + W(sin$i + sin $ 2 ) 

+ ia 2 (cos 2 $1 + cos 2 $ 2 ). (48) 

This reduced potential reveals an interesting formal anal- 
ogy to the case of an easy-axis antiferromagnet with ex- 
change constant p 2 , anisotropy a 2 , and an effective field 
of strength hd applied along the easy b axis; even though 
the actual field points along the c axis. 

Therefore the search for the minima of (48) follows the 
familiar pattern of the conventional spin-flop transition in 
an easy-axis classical antiferromagnet, with due attention 
to the fact that the relevant order parameters are now 
the staggered moments and not the actual spins. The 
results of this straightforward minimization problem can 
be simply stated by introducing the temporary notational 
abbreviations 

x=™ y = % (49) 

p p 

and are summarized in the T = phase diagram of Figure 
4 supplemented by the typical configurations within each 
phase illustrated in Figure 5. 
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hd/p 



FIG. 4. The T = phase diagram for a field h parallel 
to the c axis. The true critical boundaries are depicted by 
solid lines, the limit of local stability of Phase I is shown by a 
dashed line, and the same limit of Phase III by a dotted line. 
The WF transition in La^CuOA is described by the first-order 
I : III transition. 



A a 



A c 





a 















b 


n 2 f 







II 



11 



m i t t ni 



l I I "'2 



III 



m, I T m 



FIG. 5. Representative ground-state configurations in 
Phases I, II, and III of the T = phase diagram of Figure 
4. In all three cases the net magnetization m = |(mi + 1112) 
points along the field direction (c axis). 



At zero field (X — 0) the ground state is given by equa- 
tion (45) and exhibits twofold degeneracy. It is sufficient 
to consider the configuration defined by the upper sign 
in (45) and follow its evolution at nonvanishing field h: 



ni = -e&, 
n 2 = +e b , 



mi = \(h + d)e c , 
m 2 = 4(ft-d)e c , 



m = i(mi + m 2 ) = f fte c 



(50) 



which is depicted in entry I of Figure 5 and exhibits a net 
moment m along the c axis whose magnitude increases 
linearly with the applied field. 

Configuration I remains locally stable until the field 
crosses the boundary 



Y = y/l + X 2 -1 



(51) 



shown by a dashed line in the phase diagram of Figure 
4. However this state becomes metastable at an earlier 
stage and the true critical boundary of Phase I consists 
of two branches: 



Y = 1 - v 7 ! -X 2 , 



X < 1, 



and 



X = l, 



Y>1, 



(52) 



(53) 



which are drawn by solid lines in Figure 4 and join a third 
critical line 



X + Y 



Y<1, 



(54) 



at the "tricritical" point X = 1 = Y. 

For anisotropies below the tricritical point (Y < 1) the 
system would undergo a first-order transition at the crit- 
ical line (52) to enter Phase II characterized by a flopped 
configuration of the staggered moments ni and n 2 but 
magnetizations mi and m 2 that are both aligned along 
the c axis. With further increase of the applied field 
beyond the critical boundary (54) a second-order tran- 
sition occurs and the system enters Phase III in which 
both staggered moments are parallel to the (negative) b 
axis. 

As we will see in Section 4.3, the parameters of 
La-iCuO^ favor the value Y = a\/ p 2 ~ 2 and hence the 
relevant weak-ferromagnetic (WF) transition is the di- 
rect I : III transition that occurs at the critical line (53), 
or at a critical field h = h given by 



ho = p 2 /d. 



(55) 



The WF transition is clearly first order because the 
boundary of local stability of Phase I shown by a dashed 
line in Figure 4 extends well to the right of the true crit- 
ical boundary (53). Similarly the boundary of local sta- 
bility of Phase III shown by a dotted line extends well to 
the left of the true critical boundary. 
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One may express equation (55) in terms of the original 
variables to write 

mH = s 2 J ± , J x = 4( J x - J a ) , (56) 

where Hq is the critical field, m is the weak moment 
per Cu atom at zero field, s — 1/2 is the spin of a 
Cu 2+ ion, and J± is an effective interlayer exchange 
constant. With this identification of J± equation (56) 
coincides with the original estimate of Thio et al. |Q. 
In oxygen-doped La2CuOi+ y samples of reduced Neel 
temperature (T/v ~ 2A0K) the measured critical field is 
Hq = 5T, m = 2.1xl0~ 3 /is, and thus J± is estimated to 
be 2.5fieV. 

The description of the general features of the WF tran- 
sition is completed by quoting the explicit ground-state 
configuration in Phase III: 

ni = -e b = n 2 , 

m=mi = m 2 = -((i + /i)e c , (57) 

where the net moment m again increases linearly with h. 
The results for the net magnetization in Phases I and III, 
given by equations (50) and (57), are shown by a dashed 
line in Figure 6 which exhibits a sudden jump at the 
critical field Hq = 5T due to the first-order nature of the 
WF transition. One should stress that this calculation 
is strictly valid at T = 0; the magnetization jump is 
increasingly smoothed out with rising temperature [^|,^| . 

0.01 I ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 




magnetic field (T) 

FIG. 6. Net magnetization induced by an applied field. 
The dashed line corresponds to a field along the c axis and 
demonstrates the WF transition at the critical field Ho = 5T. 
The solid line corresponds to a field along the b axis and 
displays discontinuities at the critical fields Hi = 10T and 
H2 — 20T characteristic of the SF transition. In either case 
the net magnetization points along the direction of the ap- 
plied field. 



It is now interesting to recalculate the magnon gaps 
in the presence of the interlayer coupling and follow 
them through the WF transition. If we parametrize 
ri! and n 2 by two replicas of the angular variables (33) 
the ground-state configuration in Phase I is given by 
(9i = §,$! = -§) and (0 2 = § ,$ 2 = §)■ Small fluc- 
tuations are then studied by introducing the variables 
(0i = § = -f + <j>i) and (0 2 = f -0 a ,*2 = 

§ — 02 ) in the complete Lagrangian (43) and keeping 
terms that are at most quadratic : 

l « \(4>l + $ + + el) - p 2 (^i02 + W 
-\{ct4>l + ^4>l + 4el + c^e 2 2 ), 

cf = a 2 ±hd + p 2 ,cf =a 2 2 ±hd + h 2 + p 2 . (58) 

Standard diagonalization of this quadratic form yields 
the four magnon gaps 

n 2 1± =aj + p 2 ± ^Jh 2 d 2 + p\ 



n 2 2± =a 2 + p 2 ± ^fh 2 d 2 +p A + h 2 , (59) 

which reduce to the gaps of equation (36) at zero inter- 
layer coupling (p 2 = 0) . The zero- field limit of equation 
(59) is also interesting and leads to 

Hi- — a li ^2- = a 2, 

n 1+ = <Ja 2 + 2p 2 , n 2+ = ^a 2 + 2p 2 . (60) 

The "acoustical" gaps V and 2 _ do not depend on the 
interlayer coupling and correspond to the usual antiferro- 
magnetic (AF) modes for in-plane and out-of-plane fluc- 
tuations, respectively. The "optical" gaps 1 + and 2 + are 
sensitive to the interlayer coupling and may be said to 
correspond to exchange (E) modes, in analogy with a 
similar distinction made within a proper four-sublattice 
formulation of orthoferrites ]2l| . In the latter case, AF 
and E modes are widely separated thanks to a strong in- 
terlayer exchange interaction that is comparable to the 
intralayer one. In contrast, a close proximity of these two 
types of modes should be expected in LaiCuO^ because 
2p 2 ^a 2 . 

The decoupling of AF and E modes suggested by equa- 
tion (60) no longer holds in the presence of an applied 
field, as is apparent in equation (59). In this respect, it 
is also useful to follow the gaps beyond the WF transi- 
tion. In Phase III the ground-state configuration is given 
by 0i = § = ©2, $1 = -f = $2 and small fluctua- 
tions lead to a quadratic Lagrangian similar to (58). The 
corresponding magnon gaps are found to be 

=a 2 l +hd-2p 2 , 
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k c 



^2- 



hd + h 2 - 2p 2 



Q 2 + = a\ + hd, 



2 + 



hd+h 2 



(61) 



where the role of acoustical and optical modes is clearly 
interchanged. To be sure, equation (59) is valid for 
h < ho and equation (61) for for h > ho- The spec- 
trum exhibits a discontinuity at h — ho thanks again to 
the first-order nature of the WF transition. 

The calculation of dynamic susceptibilities can be ef- 
fected by the method briefly outlined in the concluding 
paragraph of Section 3 within the simpler 2D context. 
We do not present here explicit results in order to avoid 
further proliferation of algebraic formulas. But such re- 
sults can be provided should the need arise. 



B. Spin- flop transition 

The case of a field applied along some direction in the 
basal plane is equally interesting. In particular, when 
the field is precisely aligned with the b axis, an unusual 
spin-flop (SF) transition is reflected in magnetoresistance 
measurements 0. An important element in the cor- 
responding theoretical analysis is that the observed in- 
plane magnon gap is smaller than the out-of-plane gap 
(ax < a 2 ). 

The total bilayer potential in a field h = (0, h, 0) is 
given by 

V = /9 2 (ni • n 2 ) - hd(n lc + n 2c ) + \h 2 {n 2 b + n 2 b ) 



[«?(»?. 



L 2a) 



It 2 

■ a 2 (n lc 



4)1 



(62) 



Inspite first appearances minimization of the above po- 
tential is straightforward. The two unit vectors ni and 
ri2 are parametrized in terms of two sets of angular vari- 
ables (0i, 3>i) and (62, $2). One can then show that the 
minima of (62) are such that 



9! = e 2 = e 



(63) 



for any value of the applied field. However the azimuthal 
angles $1 and <& 2 display different behavior in three dis- 
tinct field regions separated by two critical fields: 



hi = ai 



1 



h 2 = -(a 2 



2 ' 2p 2 -a\) 



(64) 



For h < hi the ground-state configuration is illustrated 
in the first entry of Figure 7. The staggered moments 
are both contained in the be plane and cant toward the c 
axis with which they form an angle O given by 



h<h 



1 

f 




k C 



h>h- 



FIG. 7. Representative ground-state configurations in the 
three distinct field regions that characterize the SF transition. 
In all three cases the net magnetization m = |(mi + m2) 
points in the direction of the applied field (b axis). 



cos© = 



hd 



2p 2 - h 2 ' 



(65) 



Accordingly the net magnetization points along the di- 
rection of the applied field, namely 

m = — (mi + ma) = —(d + h cos 0) cos e^. (66) 

For hi < h < h 2 the staggered moments flop into the ac 
plane and form an angle with the c axis given by 



cos© 



hd 



2p 2 



while the corresponding net magnetization is 
m = mi = m.2 = —(h + c?cos©)eb. 



(67) 



(68) 



Finally, for h > h 2 , both staggered moments are aligned 
with the c axis but the net magnetization 



m = mi = m 2 = -id + h)e b 
4 



(69) 



continues to point along the field direction. This picture 
should be completed with the remark that the calculated 
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sharp SF transition at the critical fields hi and h% is 
smoothed out when the direction of the applied field de- 
parts from the b axis. 

The preceding description of the SF transition con- 
firms the theoretical analysis of Thio et al. which was 
in turn shown to be consistent with experiment fl. In 
particular, the critical fields (64) agree with those of ref- 
erence 0] if wc adopt the minimal choice a± — d and 
identify the effective interlayer exchange constant Jj_ as 
in equation (56). The net magnetization calculated from 
equations (66), (68) and (69) is depicted by a solid line 
in Figure 6 and displays characteristic discontinuities at 
the observed critical fields Hi = 10T and H 2 = 20T. A 
minor difference in the overall scale of Figure 6 with the 
corresponding result of reference |5| is apparently due to 
a slightly different choice of parameters, as discussed in 
the following subsection. 



C. Rationalized units and constants 

Our purpose here is to demonstrate how to efficiently 
use the rationalized formulas derived throughout this pa- 
per, rather than to analyze in depth the available experi- 
mental data. Such an analysis is complicated by the fact 
that actual experiments have been performed on sam- 
ples with varying oxygen doping. Pure La 2 CuO± samples 
have been available [jl 32 1 and exhibit magnetic order be- 
low the Neel temperature T/v = 325X. However the most 
complete set of magnetic measurements was obtained on 
oxygen-doped La 2 Cu04 +v with reduced Neel tempera- 
ture P|||- Hence our demonstration will be based on 
the latter measurements but could be extended to pure 
samples in a straightforward manner. 

We begin with a parameter-free theoretical prediction 
based on the fact that the first critical field in equa- 
tion (64) and the zero-field 1~ gap in equation (60) are 
both equal to a±. In physical units this equality reads 
9m^Bii\ = Qi— where g m = 2.2 is the gyromagnetic ra- 
tio and hb the Bohr magneton. Hence, if we use the 
measured critical field Hi = 10.5 ± IT, the predicted in- 
plane gap Hi- = 1-33 ± 0.l2meV is consistent with the 
measured 1.1 ± 0.3meV. For simplicity we adopt in the 
following the rounded critical field value Hi = 10T which 
leads to f^- = \.27meV . 

Now the theoretical zero-field weak moment per Cu 
atom is m = 8d/A or, in physical units, m = sgmpsSd/A 
which should be compared to a measured value 2.2 x 
10~ 3 /ib to yield 6d = 8 x 10~ 3 . This is an estimate of 
the DM anisotropy recalling that 8d = 2D /J. However 
a more convenient framework is obtained by exploiting 
the scale parameter 6 to define rationalized units such 
that d = 1, and thus S = 8 x 10 -3 , as anticipated by 
the discussion of Section 3. For the moment we restrict 
attention to the minimal model for which ai = d = 1. 
Hence the theoretical critical field hi = a\ = 1 sets the 



rationalized field unit equal to the measured Hi = 10T, 
and the theoretical in-plane gap Qi- = a\ = 1 sets the 
unit of frequency at 1.27meV '. Then the measured criti- 
cal field for the WF transition Hq — 5T is translated into 
ho = 1/2 rationalized units and thus equation (55) reads 
ho = p 2 /d = p 2 = 1/2 which provides a rationalized esti- 
mate of the interlayer coupling. Finally we consider the 
theoretical critical field h 2 of equation (64) which may 
now be applied with a\ = d = 1 and 2p 2 = 1 to yield 
h 2 = a 2 , in rationalized units or H 2 = 10a 2 T in physical 
units. Comparing this prediction to the measured critical 
field H 2 = 20T we find that a 2 = y/2. 

To summarize, all theoretical formulas may be applied 
with rationalized parameters 



5 = 8 x 10~ 3 ; 



a i 



1. 



a 2 



V2, 



(70) 



supplemented by the stipulation that the physical unit 
for frequency be l^lmeV, for field 10T, and for magne- 
tization sg-mP-B = 1.1/XB. 

The predicted values for the acoustical gaps f^- = 
1.27<zi = 1.27meV and £1 2 - = 1.27a 2 = 1.8mey are 
marginally consistent with the observed 1.1 ± 0.3meV 
and 2.5 ± 0.5meV. In fact, full consistency would be 
restored if we had included error bars in our analysis j|] 
instead of the conveniently rounded input values for the 
critical fields and the magnetization actually used in our 
demonstration. Furthermore the optical gaps of equation 
(60) are now predicted to be 0] + = 1.27 '\Ja\ + 2p 2 — 
l.8meV and fl 2 + = 1.27 y/a% + 2p 2 = 2.2meV. We do 
not know of an experimental determination of the optical 
gaps. Thus we merely note the predicted close proximity 
of acoustical and optical gaps, as anticipated in Section 
4.1. 

Since the unit of frequency is equal to 2s5J = 5 J — 
1.27meV, the exchange constant is predicted to be J = 
1.27/i5 = l58meV. Curiously, this classical value of the 
exchange constant is the same with the one obtained in 
Section 3 in relation to the spinwave velocity 850meVA 
observed on pure (T/v = 325if ) samples. But the spin- 
wave velocity on oxygen-doped samples is typically lower 
(~ 700meV r A) and thus the currently predicted classi- 
cal J is somewhat uncomfortably high. One would think 
that such a discrepancy can be averted by resorting to 
a nonminimal model with ai ^ d. Interestingly, a more 
comfortable value of J can thus be obtained but only at 
the expense of further deterioration (lowering) of the pre- 
diction for the out-of-plane acoustical gap fl 2 - — 1.8meV 
discussed in the preceding paragraph; and vice versa. 
Putting it differently, current data do not indicate de- 
parture from the minimal model ai = d. 

The overall consistency is satisfactory in view of the 
following limiting circumstances: 

(a) The theoretical model is strictly valid for pure 
La 2 CuOi and its application to oxygen-doped samples 
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should be viewed with caution. Unfortunately, a com- 
plete set of magnetic measurements is not available on 
pure samples probably because the observed magnon 
gaps are nearly doubled H and thus the critical fields re- 
quired to induce the SF transition become prohibitively 
intense. 

(b) The theoretical calculation is strictly valid at T = 
0. Although the measured magnetization and critical 
fields have been extrapolated to T = 0, other quan- 
tities such as the magnon gaps are typically determined 
at relatively high temperatures (T ~ 100K). 

(c) The assumption of strictly localized spins is always 
in question, and the further use of a classical description 
is marginal for this low spin value (s = 1/2) and low ef- 
fective lattice dimension (D = 2). Since different physical 
quantities are renormalized differently by quantum fluc- 
tuations, including calculated quantum-renormalization 
factors into the classical predictions can be confusing [R2J . 



restoring force supplied by the antiferromagnetic inter- 
layer coupling. When the field exceeds H stripes are 
rendered unstable and begin to expand steadily in both 
directions, thus providing a detailed mechanism for the 
observed first-order WF transition. 

The preceding qualitative picture may be put on a firm 
quantitative basis using the derived continuum field the- 
ory, as we hope to demonstrate in a future publication. 
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V. CONCLUSION 

While the phenomenological picture derived by Thio et 
al. f|-§| is confirmed by the present analysis, some new 
elements have emerged that may deserve closer attention: 

The structure of the magnon gaps is more involved 
than normally assumed because of the underlying four- 
sublattice magnetic ground state. The calculated acous- 
tical and optical gaps are not widely separated and hy- 
bridization takes place in the presence of external fields. 
Therefore study of the field-dependence of the magnon 
gaps may lead to further tests of the derived picture. 

The description of the isotropic 2D antiferromagnet 
in terms of a relativistic nonlinear a model has already 
provided interesting results |l||Jl6| but the presence of 
anisotropics and an interlayer coupling are clearly impor- 
tant for a more detailed understanding of the magnetic 
structure of La-iCuO^. Suffice it to say that the existence 
of a finite Neel temperature (Tjv = 325K) is precisely due 
to such perturbations. Hence it is conceivable that the 
field theoretical framework discussed here may help to 
address some of the remaining questions J32| . 

The covert nature of weak ferromagnetism in LaiCuO^ 
makes it difficult to directly observe macroscopic mag- 
netic domains. Nevertheless, even on structurally pure 
samples, cooling below Tjv should produce numerous 
magnetic domains and antidomains separated by domain 
walls that are invisible because the average magnetiza- 
tion vanishes at zero external field. When a field is ap- 
plied in a direction perpendicular to the CuOi planes 
domain walls evolve into magnetic stripes that exhibit 
enhanced magnetization over a region of finite width and 
could thus become visible. Static magnetic stripes are 
stable for field strengths smaller than the critical value 
Hq required to induce the WF transition, thanks to the 



APPENDIX A: NNN IN-PLANE INTERACTIONS 

Here we consider the modifications of the 2D Hamilto- 
nian of Section 2 that result from spin interactions along 
the diagonals of the Cu plaquettes. Symmetry precludes 
the existence of antisymmetric DM anisotropies on such 
bonds and thus all nnn contributions to the Hamiltonian 
may be cast in the form 



«fcZ> i,j 



(Al) 



where <C kl 3> denotes a nnn bond and the symmetric 
matrices Gu = (G^) are not assumed to be traceless. 
The possible value of Gki are again restricted by sym- 
metry as shown in Figure 8. Thus A spins interact with 
their A neighbors through the exchange matrices 



Gi = 



G 



G 2 



laa 











Gibb 


Guc 





Gibe 


Glee 


laa 











Gibb 


Gibe 





Gibe 


Glee 



(A2) 



along the a and b direction, respectively, whereas B spins 
interact with exchange matrices 



G'i = 



G' 



Glaa 











Gibb 


— Gibe 





— Gibe 


Glee 


Glaa 











Gibb 


— Gibe 





— Gibe 


Glee 



(A3) 
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which are related to G\ and Gi. A corollary of these 
symmetry relations is that the isotropic components of 
nnn exchange couplings, given by the traces of the above 
matrices, are characterized by two exchange constants 
which are generally different along the a and b directions 
but the same for AA and BB bonds. 




9 = C -j-j -, G = .g(n c e fc + n fc e c ), (A4) 

where we have also defined a vector G that depends on 
the staggered moment but plays a role similar to the DM 
vector d = de a . Then equations (30) generalize to 

m=-[nx(h + d + nxG — nxh)], 

1 3 

L = -n 2 + -g(n 2 b - n 2 c )n a + h • (n x h), 

V = (h x d) • n + g [h c rib + /i&n c — 2(n • h)nbn c ] 

+ ^(n • h) 2 + l -{a\n\ + a 2 2 n 2 c ) + 2g 2 n 2 b n 2 c , (A5) 

which were actually used throughout our analysis. But 
the results presented in the main text were restricted to 
g = mainly because the qualitative picture remains 
intact and current experimental data are not sufficiently 
accurate or detailed to discern a nonvanishing g. 




'a,p-l A a+l,p B a+l,p 



FIG. 8. Illustration of the distribution of the matrices G 
and G' given in the Appendix. These symmetric matrices 
describe nnn exchange interactions along the diagonals of the 
Cu plaquettes which are parallel to the a and b axes. The 
remaining conventions are those of Figure 2. 

The corresponding modifications of the effective low- 
frequency dynamics may be briefly summarized as fol- 
lows. The traces of the above matrices introduce an 
overall additive renormalization of the isotropic exchange 
constant J. Similarly the contributions from the diago- 
nal anisotropies submerge with the anisotropy constants 
a\ and a 2 already discussed in the main text. The only 
new parameter is then introduced by the off-diagonal 
anisotropy, namely 
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